Onsager-Manning-Oosawa condensation phenomenon and the effect of salt 



Emmanuel Trizac 

CTBP, UC San Diego, La Jolla (USA) and LPTMS Universite Pans XI, 91405 Orsay (France) 



in 
o 
o 

O 

Q 



Gabriel Tellez 

Departamento de Fisica, Universidad de Los Andes, A. A. 4976, Bogota, Colombia 

Making use of results pertaining to Painleve III type equations, we revisit the celebrated Onsager- 
Manning-Oosawa condensation phenomenon for charged stiff linear polymers, in the mean-field 
approximation with salt. We obtain analytically the associated critical line charge density, and 
show that it is severely affected by finite salt effects, whereas previous results focused on the no 
salt limit. In addition, we obtain explicit expressions for the condensate thickness and the electric 
potential. The case of asymmetric electrolytes is also briefly addressed. 
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Whereas scaling approaches have proven useful to de- 
scribe neutral polymers, our understanding of polyelec- 
trolytes solutions is quite rudimentary, due to the long 
range character of Coulombic interactions 0,0. Stiff lin- 
ear polyelectrolytes that are free of the coupling between 
chain conformation and small ions degrees of freedom 
consequently provide ideal systems for a comprehensive 
comparison between theory and experiments. There in- 
deed exists a large variety of such rod-like polyions, rang- 
ing from synthetic polymers (e.g those based on poly(p- 
phenylene) backbones to biological molecules (DNA, 
actin filaments, micro-tubules, some viruses etc). On dis- 
tances smaller than their (large) persistence length, these 
objects behave as charged cylinders with an associated 
logarithmic electrostatic potential that may be strong 
enough to bind oppositely charged microions (counteri- 
ons). This was first realized by Onsager and analyzed 
by Manning £| and Oosawa p| : in the limit of van- 
ishing polymer radius, the corresponding phenomenon of 
counterion condensation is triggered when the so-called 
Manning parameter £ = £b/£ is larger than unity (see 
e -g- UL Q f° r the essence of the argument). Here, i~ x 
is the backbone line charge in units of the elementary 
charge e and Ib = e 2 /(ekT) denotes the Bjcrrum length 
with e the dielectric constant of the solvent and kT the 
thermal energy. 

Counterion condensation, which affects a gamut of 
static and dynamical properties, is central to our view 
of highly charged polymers and an active field of re- 
search [i 0, H H OB, (HI- Interestingly, the influential 
arguments of Onsager, Manning and Oosawa found con- 
firmation and a firm basis in numerical and analytical 
studies of the mean-field non-linear Poisson-Boltzmann 
(PB) equation where the densities of microions with va- 
lency z are related to the local mean electrostatic poten- 
tial ip by n z cx exp(—zetp/kT). The basic features are 
already present in the exact solution of PB theory for a 
charged cylinder in a concentric cylindrical Wigner-Seitz 
cell without salt which can be cast as a restricted 
version of a partial differential equation first studied and 
solved by Liouville in 1853 [l^, as pointed out in [T3 |. 



For an isolated cylinder (i.e. in the limit of a diverging 
Wigner-Seitz cell radius) with a small charge £ < 1, no 
counterions are bound to the cylinder and the potential 
distribution reduces to the bare logarithmic form. On 
the other hand, for £ > 1, a finite fraction of ions are 
"bound" to the surface |l5j. 

The physically relevant situation where salt is present 
(i.e. with both co- and counter-ions, that may come from 
the dissociation of an added electrolyte and/or from the 
solvent itself) has resisted analytical understanding much 
longer. Ramanathan has showed that for £ > 1 and dis- 
tances much larger than the polyion radius a the ionic 
atmosphere was the same as that due to a cylinder with 
charge parameter £ = 1 [lq . This result holds asymptot- 
ically for na — > |l6t fl7| where k _1 is the Debye length 
pi • This condensation effect is the counterpart of the 
aforementioned one and was placed within the mathe- 
matical framework of isomonodromy theory in Ref. [ts| 
and Painleve III equations in Refs. [3 [l! H(j. This 
drew an exact correspondence with Ising model corre- 
lators but more importantly allowed a) to compute the 
exact far-field behaviour [Ig> b) to show rigorously the 
critical nature of the value £ = 1 when na — > |^.ll8j. and 
c) to obtain some analytical results for the potential dis- 
tribution and ionic densities Q . The common framework 
to the previous approaches is PB theory [2l| to which we 
will restrict ourselves in the subsequent analysis. Most of 
existing results hold for na — > only, and it appears that 
this (singular) limit is approached logarithmically slowly. 
From a practical point of view, finite na corrections are 
therefore crucial and can never be discarded. 

The purpose of the present work is to analyze the 
fate of the counterion condensation phenomenon at fi- 
nite salt concentration {na ^ 0). We shall show that for 
na < 1, the transition is smoothed but remains, and that 
the associated critical charge parameter is salt-dependent 
and smaller than the usual Manning threshold £ c = 1. 
The condensate structure will be resolved without the 
need to invoke any matching procedure, and analytically 
tractable results will be obtained for the electric potential 
below and above the critical value £ c . Extensive use will 
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be made of results pertaining to the theory of Painleve 
III equations [2(|, that mathematical difficulty have - 
to the best of our knowledge- prevented to find their 
way towards the physicists' community j^. The case of 
asymmetric f:2 and 2:1 electrolytes -that are physically 
not equivalent [HI - w ^ a l so be briefly addressed. Most 
of our results hold for na < 1 and by comparison with 
numerical results can be shown to be very accurate in 
this parameter range. 

We consider Poisson's equation in cylindrical coor- 
dinates for the dimensionless electric potential [(f) — 
e^/(kT)} : 
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for the situations (z + ,z_) = (1,— 1), (1,-2) or (2, — f), 
other asymmetries unfortunately resisting the analysis. 
We demand that rdip/dr = — 2£ for r = a where the 
dimensionless bare charge of the polyion, £, is assumed 
positive without loss of generality. The second boundary 
condition (0 — *■ when r — > oo) ensures that at large 
distances, tf> obeys the linearized equation V 2 </> = n 2 (j) 
and therefore takes the form 
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where Kq and K\ denote the modified Bessel functions of 
the second kind, of order and f respectively. The pref- 
actor defining the a priori unknown effective charge £ c ff 
reduces by definition to the bare charge £ for low values of 
£ (linear regime) but non-linear effects significantly affect 
£ c ff. In particular, a consequence of countcrion conden- 
sation is that for na — > 0, the former quantity no longer 
depends on £ provided £ > 1 [H 0, E| • We first focus 
onto 1:1 electrolytes for which the electr ic p otential <j>u 
crucially depends on a parameter A 0, that in the 
present context is related to the effective charge through 



£ ff = 2 na K\(k<x) A 



(3) 



so that at large distances, <f> ~ -IXKo(Kr). For A < 1/tt 
(which corresponds to £ < £ c where £ c is a critical charge 
to be defined below), theorem 3 of Ref. [2(| implies the 
short distance behaviour 
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r being the Euler function. On the other hand, for A > 
1/tt, there is a qualitative change of behaviour for (f>n 
that then reads 
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where denotes the argument of T(ifi) J20|. For our 
purposes, a very accurate approximation -considered in 
the remainder- is given by the small \x Taylor expansion 
<f(/i) ~ -7r/2 - 7/i, with 7 ~ 0.5772... the Euler con- 
stant. Both quantities a (defined only for A < 1/tt) and 
/i (defined only for A > 1/tt) are positive and related to 
the bare charge by the boundary condition r<p'(r) = — 2£ 
at r = a : expressions (@J and |J?jJ yield 
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(£-l)tan[2/xlog(/to/8) + 2/x7] = 
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where it is understood that /i is the smallest positive root 
of |0. As such, a and ^ characterize the short distance 
features, and the difficulty amounts to connecting these 
parameters with the far-field quantity A. It may be shown 

that miu 
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for A < TT' 1 ( or £ < £ c ) (9) 
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A = - cosh(7r^) for A > 7T ( or £ > £ c ) (10) 

TT 

The leading order term in cf>n is therefore — 2crlog(Kr) 
for A < 7T _1 which corresponds to the bare potential of 
a cylinder with line charge a. For A > tt^ 1 the dom- 
inant small r behaviour reads — 21og(«;r), up to an r- 
independent term varying with charge and salt content 

M. 




FIG. 1: Contour lines for a and fi (and hence for A) in the 
charge-salt plane, as deduced from Eqs JJJ and JHJ. The thick 
dashed curve shows the locus of critical charges £ c and sepa- 
rates the iso-<r curves for £ < £ c (from left to right, the curves 
correspond to a — 0.2, 0.4, 0.6 and 0.8) from the iso-n curves 
shown -on the right hand side- for £ > £ c (from bottom to 
top, (j, = 0.08, 0.1, 0.13 and 0.2). 

The change of behaviour displayed by Eqs. JI) and 
© is the fingerprint of counterion condensation. The 
exponent a fulfills < a < 1 and for the critical value 
A = 1/7T, we have a — 1 while /i vanishes. The corre- 
sponding critical value of £ follows either from @ taking 
the limit a — > 1~ or enforcing fi — ► + in Eq. (JSJ. These 
two routes yield exactly the same critical charge (which 
illustrates the consistency of the underlying expressions 
for the potential): 

ic = 1 + logM+7-log8- 
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In the limit Ka — > 0, the celebrated Manning threshold 
£ c = 1 is recovered but the correction embedded in 
is significant: at na — 10~ 3 , £ c is shifted to 0.881 and at 
Ka = 0.1, we get £ c ~ 0.737. Expression (fill) has been 
derived omitting corrections of order (kt) 2 in and of 
order (ftr) 4 in ©. One may show however that the next 
correction to (|llf> is of order (na) 4 (log Ka) 2 , and is there- 
fore irrelevant from a practical point of view whenever 
na < 1 [2!| . Figure may be considered as illustrat- 
ing a "law of corresponding states" and shows the values 
of bare charge and salt concentration (or equivalcntly 
polyion radius) leading to the same electrostatic poten- 
tial: if a or fi is fixed, the connection formulae © and 
(fT7Tf> indeed ensure that A is also fixed, so that moving 
along the contour lines of Fig leaves the full function 
4>(nr) unaffected. The complementary information of the 
£ dependence of a and fi at fixed salt concentration is 
shown in Figure |2 It appears that except in the vicinity 
of £ c , cr is very close to £ (with systematically a > £). 
Close to the transition threshold, one has a — 1 and fi 
scaling like |£ — Cd 1 ^ 2 - For highly charged polyions, the 
argument of the tangent function in (JSJ) has to be close 
to — 7r, from which we obtain 

-7I-/2 

M ~ log( Ka )+ 7 -log8- (£- (1 ' 

The inset of Fig. [3 assesses the quality of this approxi- 
mation which improves when £ increases but fails quite 
severely when £< 1 + 0(1/ log Ka). 
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FIG. 2: a versus bare charge. In the main graph and in the 
inset, 3 values of Ka have been chosen. From left to right, 
Ka = 10" 2 ,10" 3 and 10" 6 . For Ka -> 0, a -> £ (see the 
dashed first bissectrix) while /j, — > 0. The associated critical 
charges are indicated by the arrows (also in the inset). The 
inset shows fi versus £ and the dashed curves correspond to 
approximation 1121 . The vertical line indicates f = 1. 

We have tested the accuracy of expressions (@J /Q)/ © 
below the critical charge, and Eqs. ©/©/CSl f° r 
£ > £ c against numerical solutions of the full non-linear 
PB equation lJT|l. These formulae turn out exception- 
ally accurate for small r, while the far-field is conversely 
equally well given by where £ c ff is determined by J3J). 
More interestingly, for kt < 1/2, <fiu deduced from Q is 
extremely accurate for all possible charges and Ka < 0.1 
(curves not shown). Similarly, the far-field expression 
@ provides very good accuracy as soon as kt > 1/2. It 



therefore appears that <j)n is analytically known at all 
distances ; in the crossover region kt ~ 0.5 where no 
limiting expression is supposed to be reliable, the worst 
relative accuracy is observed but is nevertheless always 
below 2%. We emphasize here that £ e ff significantly dif- 
fers from the critical charge £ c . In other words, even 
treating correctly the condensation phenomenon, non lin- 
ear effects are still present and affect the diffuse cloud of 
remaining "uncondensed" ions, that cannot be treated 
within linearized PB or Debye-Hiickel theory. This im- 
portant aspect has been overlooked so far (see e.g. 0) 
and is present even in the simplest limit Ka — > 0, where 
£ c ff —> 2/-7T irrespective of £ (provided £ > £ c ), whereas 
£ c = 1. The resulting repulsive force of interaction be- 
tween two highly charged rods is therefore overestimated 
by a factor T — (tt/2) 2 ~ 2.5 in the traditional Manning- 
Oosawa picture. An experimental check of this prediction 
has to fulfill the requirement of very low salt, since at say 
Ka = 10 -2 , the above ratio is closer to unity [T ~ 1.4 as 
can be checked from Eq. ©]. 

The previous results allow to discuss in a rigorous 
manner several aspects of polyelectrolyte physics inves- 
tigated in the literature. Ofparticular interest is the 
so-called Manning radius Rm P.ITol|. introduced to quan- 
tify the condensate thickness when £ > £ c . In the 
Wigner-Seitz cell approach, the integrated line charge 
per Bjerrum length £b [which, from Gauss' law reads 
q(r) = —rcj)'(r)/2] has an inflexion point when plotted as 
a function of logr at the distance Rm where q(Rn) = 1 
[ji|2ij. It is remarkable that the q(r) following from © 
shows the very same feature, which is ascribable to a sim- 
ilar functional form, although at finite salt concentration 
q(a) = £ = 1 no longer provides the critical value of the 
charge. The corresponding Manning radius reads: 

kR m = 8exp[-7r/(4/x)- 7 .] (13) 

When £ > 1 + C(1/ log Ka) , one may use the approxima- 
tion (I12|) and therefore 

kR m ~ 2V2^ exp (-1- 2 * ) . (14) 

The latter formula is the central result of a recent work 
[To| where it was obtained by a clever but approximate 
matching procedure. It appears to be incorrect when 
H deviates from the large £ expansion l|12l) (see the in- 
set of Fig. [2] where it is seen that the difference be- 
tween the dashed and continuous curves may be large). 
On the other hand, Eq. Ijl3(l is always found to be 
very accurate compared to the numerical solutions of PB 
equation [2j|. We note that choosing a different defi- 
nition for the condensate thickness R* (e.g. demand- 
ing that q{R*) — £ c instead of unity, leads for large 
bare charges to the scaling relation R* oc a a K a ~ 1 with 
now a — (arctan7r)/7r ~ 0.402, smaller than the value 
a = 1/2 appearing in i|14|) . The convenient inflection 
point criterion alluded to earlier would however be lost 
following such a route. 
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We now turn to the case of asymmetric 1:2 and 2:1 
electrolytes [23|. Recent results obtained for a class of 
solutions to the cylindrical Toda equations [23] allow to 
extend the above analysis to such situations. So far, the 
1:2 case only has been studied, again only in the limit 
na — » 8j. The details will be provided elsewhere pfij : 
we concentrate here on the critical charges 

1 



^ 2 + 21ogM + 2C 1:2 (15) 

c A = i + , ! „, , (is) 

log(Ka) +C 2:1 v ; 

with C 1:2 =7 - (log2)/3 - 3(log3)/2 ~ -1.301(17) 
and C 2:1 = 7- log 2 - 3 (log 3) /2 ~ -1.763. (18) 

The associated Manning radii, defined from g(i?Af) = 
1 (2:1 case) or q(Rm) — 1/2 (1:2 case) -which again 
ensures the existence of an inflection point criterion- read 



KR\f ~ 2~ 1/3 3 3/4 V2^a~exp 



kR 2:1 



exp 
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(19) 
(20) 



These expressions have the same status as (|14|) . i.e. hold 
at high enough £. It appears that the Manning 2:1 ra- 
dius is inflated a factor 3 3 / 4 /2 ~ 1.14 compared to its 
1:1 counterpart, irrespective of salt content and charge 
(but beyond the condensation threshold). This quanti- 
fies the intuitive picture of a swollen double-layer due 
to the presence of divalent coions (expelled further away 
than monovalent ones) , and conversely of a shrunk cloud 
(by a factor 2~ 4 / 3 3 3 / 4 ~ 0.90 when £ is large enough) 
in the 1:2 case due to the more efficient screening with 
divalent counterions. 

To conclude, exploiting an important body of mathe- 
matical work in the field of stiff polyelectrolytes allows 
to systematically address finite salt effects, that are cru- 
cial even under experimentally low salt conditions. The 
present work sheds new light on the condensate structure 
and further connects short scale features with the long 
range behaviour of the electric potential. Among the con- 
sequences of experimental relevance following from our 
analysis, we mention the large q behaviour of the partial 
counterion/rod S- and coion/rod S+ partial structure 
factors, expected to scale like q- 2±2a for £ < £ c , with 
S- ex 1/logg at £ = £ c . Such effects are missed fol- 
lowing the classical Manning-Oosawa arguments, which 
generally result in an underestimation of screening. 
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